Geometrical Frustration in Two Dimensions: Idealizations and Realizations of a Hard 

Disc Fluid in Negative Curvature 
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We examine a simple hard disc fluid with no long range interactions on the two dimensional space 
of constant negative Gaussian curvature, the hyperbolic plane. This geometry provides a natural 
mechanism by which global crystalline order is frustrated, allowing us to construct a tractable model 
of disordered monodisperse hard discs. We extend free area theory and the virial expansion to this 
regime, deriving the equation of state for the system, and compare its predictions with simulation 
near an isostatic packing in the curved space. Additionally, we investigate packing and dynamics on 
triply periodic, negatively curved surfaces with an eye toward real biological and polymeric systems. 



I. INTRODUCTION 

The true nature of glasses and the glass transition has 
long remained inscrutable. Glasses belong to a wide class 
of disordered materials and systems which include gran- 
ular materials like sand [l| or wheat grains [2I] , frustrated 
systems like spin glasses j^], and soft systems like foams 
Q and emulsions Unfortunately, the analogies af- 
forded by this classification have done little to lift the 
veil of mystery: from misunderstandings about apparent 
"flow" displayed by stained glass windows from the mid- 
dle ages to debate over the meaning of randomness 0] , 
few aspects of the problem have escaped some measure 
of controversy. The prominent role of geometrical frus- 
tration in these systems, however, has long been agreed 
upon and it is here that we turn our attention and 
expand upon our prior work Q . 

We begin by considering an idealized colloidal glass, 
made of monodisperse spheres in three dimensions with 
no long-range interactions and a purely hard-core, short- 
range potential. The entropically favored, high-density 
fee lattice (there is a virtual tie between fee and hep 
lattices, and lattices with the stacking faults that relate 
them) for the system is well understood and highly stable, 
but, nevertheless the system can exhibit a glass transition 
to a dynamically arrested, amorphous state at a packing 
fraction significantly lower than that of the crystal. One 
way to understand this is to recognize that there is a 
competition at play between local and global properties 
of the spheres' arrangements. On the one hand, the best 
the colloids can hope to do, entropically speaking, is to 
arrange themselves into one of the close-packed, global 
crystal configurations. On the other hand, these crystal 
configurations are not the best local arrangements - a 
tetrahedral packing is locally preferable, but cannot be 
maintained beyond small distances without the introduc- 
tion of defects. 

In spite of this qualitative understanding of the impor- 
tance of geometrical frustration, the full colloidal glass 
system remains difficult to tame, particularly near the 
transition [IQj . In principle, one could imagine construct- 
ing a simpler system that shares enough of the qualities 
of the colloidal glass that a full analysis yields insights 



into the more complicated case. A time-honored tradi- 
tion along these lines is the study of the same system 
in fewer dimensions [llj . In this case, we would have 
monodisperse, hard discs with no long-range interaction. 
Notice, however, that we have lost geometrical frustra- 
tion! The entropically-favored, global, hexagonal crystal 
is composed precisely of the locally optimal units, trian- 
gular configurations of discs. Indeed, two-dimensional, 
monodisperse discs do not exhibit a glass transition [lo[ . 

How then, can we use a simplified system to say any- 
thing about the genuine article? We must find a way to 
modify the two dimensional case so as to recover some 
measure of geometrical frustration. One method that has 
been well studied is the use of bi- or polydisperse discs 
(To| . frequently with disc radii in a ratio of 1.4 : 1. The 
use of two or more disc sizes successfully reintroduces 
geometrical frustration and recovers the glass transition. 
The drawbacks to this tack are the possibility of a phase 
transition into locally monodisperse, crystalline regions, 
and, the explicit introduction of additional particle types 
and the arbitrary nature of the size distribution, breaking 
the clean analogy with the original colloidal glass. 

So, is there a way we can keep the two-dimensional sys- 
tem monodisperse and yet still manage to frustrate the 
local and global packings? Indeed, in the seminal work 
of Kleman, Nelson, and others [1, [l^ , background cur- 
vature is introduced to the system, achieving the desired 
frustration. To picture how this might happen, imagine 
that our two-dimensional space is a sphere, and that the 
hard discs are embedded in it as spherical caps. Notice 
that there is no longer enough room around each spheri- 
cal cap to fit six additional caps, as there would be in flat 
space. More precisely, the number of discs (or spheres or 
hyperspheres) that can simultaneously touch a central 
one, the kissing number, nkiss, varies continuously as a 
function of the curvature, K. In two dimensions, this 
relationship is given by 13]: 
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In order for local triangular packings of discs to be com- 
mensurate with the best global packing, nkiss must be 
an integer. This condition is satisfled only for flat space 
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{K = — > n-kiss = 6), and for isolated non- vanishing 
curvatures. 

Historically, curved backgrounds in this context have 
been used to study crystalline states on curved surfaces 
[Tst . and the ordered defect patterns that emerge as 
a result, or to imagine the ideal, unfrustrated, ordered 
configurations in three dimensions with the hope of un- 
derstanding how "decurving" might lead to the physical, 
disordered states we observe in glassy systems a la the 
cholesteric blue phase [1, In contrast to these ap- 
proaches, we choose to study the disordered fluid regime 
of hard-core systems in uniform curvature. Positive cur- 
vature may, at first blush, seem more suitable due to its 
intuitive appeal and the fact that a sphere may be em- 
bedded in three dimensional Euclidean space (K^) with- 
out distortion. However, since the surface area of the 
sphere is finite, and smaller particles on the sphere's sur- 
face probe ever lower curvatures (see Appendix A) , there 
is no way to extract a sensible thermodynamic limit at a 
given curvature. 

Instead, we turn to the space of constant negative 
curvature, the hyperbolic plane, H^, where the ther- 
modynamic limit is not impossible, but remains subtle 
nonetheless. While the sphere suffers from not having 
enough room for a thermodynamic limit, on there may 
actually be too much room. To see why this might be so, 
consider a circle of radius r on H^. Its area is A{r) — 
-2nK-^ [cosh(V^r) - l] « Trr^ - f^Kr^ for smaU r, 
consistent with the known K = result. Meanwhile, the 
circle's circumference is C(r) = 27r sinh(V —Kr) /V —K^ 
which grows as fast as A{r) for large r. Thus, the area 
within an arbitrary distance of the boundary of a circle, 
or any other fixed shape, is always a finite fraction of 
the entire region's area, no matter how large the region 
becomes. Since the traditional thermodynamic limit de- 
pends on the effect of the boundary vanishing at infinite 
size, we must tread carefully. Because our system has 
interactions that are purely local, the ideal gas provides 
a useful test case for managing the issue. 

The textbook approach to the ideal gas law begins with 
particle-in-a-box eigenenergies and takes the large area 
limit to calculate the partition function for a single par- 
ticle, Zi = V/{\tY where At = ^ 2t: /{mkBT) is the 
thermal wavelength. We might expect that the partition 
function, even on a curved manifold, would be propor- 
tional to the proper hypervolume of the box. Indeed, 
since any partition sum is dominated by the infinite num- 
ber of high energy, short wavelength terms, one would 
think that these modes are insensitive to the curvature 
as soon as their de Broglie wavelength is shorter than 
the radius of curvature. For concreteness, let's consider 
the partition function for an ideal gas on a hemisphere 
where the wavefunctions arc forced to vanish on the equa- 
tor; the energies are the eigenvalues of the Laplacian on 
the hemisphere with Dirichlet boundary conditions. On 
the whole sphere, the spherical harmonics y£„i(r2) are the 
complete set of eigenfunctions with eigenvalue £(^-1-1) for 
m = —l^ . . .1. The eigenfunctions on the hemisphere will 



be a subset of the spherical harmonics, with the bound- 
ary condition enforcing the constraint that |m| is odd 
and so the degeneracy for the l{l -I- 1) eigenvalue is I. 
Thus the partition function will be: 

where R is the radius of the spere. We can re-express this 
partition function as an expansion in At jR by employing 
the Euler-Maclaurin formula [15]: 

E/w - tmd. + l^^M (3) 

where &2fc are the even Bernoulli numbers. Substituting 
/(n) = nexp [ATn(n-|- l)/47ri?^], i = and j — oo gives 
the desired expansion: 

Notice that we get terms proportional to both the area 
and the perimeter, but also a constant piece that one 
may speculate is related to the topology of the mani- 
fold. Similarly, we could consider the partition function 
with Neumann boundary conditions along the equator, 
i.e. vanishing derivative. In this case the constraint is 
that £-f |to| is even and the appropriate sum to consider 
is 

^E(^+l)exp|-^A^| 

- ^+2A^ + 6+^('-/^) 

Again we see an area and perimeter term, the latter with 
the opposite sign, and the same constant piece. We note 
that the sum of Zi = -I- Z^ is the partition function 
for the whole sphere and we see that the perimeter terms 
cancel leaving us with 

A-kR^ 1 

Z, = ^ + - + 0{\t/R) (6) 

again, the area of the manifold plus an extra constant. 

To probe more exotic surfaces, one could attempt this 
kind of calculation on a case- by-case basis [l3|, but for- 
tunately, the spectrum of the Laplacian on a general 
Riemann surface is well characterized. Originally, Weyl 
studied the asymptotic behavior of the large eigenvalues 
Kac, in asking, "Can one hear the shape of a drum?" 
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found our constant ^ and, finally, McKean and 
Singer [l^l developed the "Weyl expansion" for the par- 
tition function of a single particle on a two-dimensional 
domain D: 
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where the sign of the first correction depends on Neu- 
mann or Dirichlet (-) boundary conditions and Kg is 
the geodesic curvature of the boundary considered rela- 
tive to its outward-directed normal [2l| . Our calculations 
are, unsurprisingly, in complete agreement with their the- 
orem which, interestingly, works in any dimension. Re- 
turning to the question at hand, namely the hyperbolic 
plane, we can see that the circumference directly enters 
the partition function at uncomfortably low order. What 
is to be done? Periodic boundary conditions provide an 
out - by eliminating the boundaries and considering suf- 
ficiently large periodic boxes we may recover some sem- 
blance of a thermodynamic limit. We hasten to point out 
that even this procedure is not free of its own subtleties 
- due to the inherent length scale provided by the cur- 
vature (Appendix A), a simple re-scaling of the periodic 
box is not allowed. In particular, by symmetry, the edges 
of the periodic box must have vanishing geodesic curva- 
ture. Thus the Gauss-Bonnet Theorem which relates the 
background curvature K, the geodesic curvature Kg and 
the jump angles at any vertices, ai, 



/ / X = 27r - V a» - / I 
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imposes a relation between the area contained within 
the box and the polygon's angles. At the same time, to 
preserve homogeneity after the periodic boundaries have 
been applied, any identified vertices must end up with 27r 
of angular space around them, thus fixing the box's area. 
Topological considerations force the number of sides to 
be a multiple of 4 [1^, and we are left with the familiar 
square cell in flat space, along with a specific-size oc- 
tagon, dodecagon, and so forth on H^. Luckily, the area 
of these p-gons increases with p , and we are reassured 
that sufficiently large periodic boxes are achievable. 

Recently, we took the first several steps in a compre- 
hensive treatment of just such a system of hard discs on 
0. Our examination included a study of ordered 
disc packings leading to the application of free-area the- 
ory, a calculation of the first several virial coefficients 
of the system, and a molecular dynamics realization of 
the fluid in the smallest possible periodic region. In this 
paper, we extend and detail these results by consider- 
ing a constant-curvature recasting of the free-area theory 
and calculating the virial coefficients to higher precision 
while applying Pade approximants to get a qualitative 
look at the behavior of the divergence in (/) as a function 
of the curvature. We expand our simulation to a larger 



polygon with more sides, allowing for better statistics 
through more simulated discs, and use it to probe the 
zero curvature limit from below. In addition, we study 
entropically governed packings on surfaces with average 
negative curvature in M? , such as the triply periodic min- 
imal surfaces found in diblock copolymers, lipids, or even 
the mitochondria of Chaos carolinensis [23j , among other 
biological examples. 

This paper is organized as follows. In section II, we re- 
view the theory of regular tesselations in two dimensions 
of arbitrary curvature, and extend this theory to packings 
of hard discs with crystalline order. Section III develops 
a free area theory from these ordered packings on curved 
space. In section IV, we calculate the first five terms 
in the virial expansion as a function of curvature, and 
comment on the analogy between high-dimensional fiat 
space and the highly curved plane. Section V presents 
a molecular dynamics model of a hard disc fiuid on 
and compares the simulation results with the predictions 
of free area theory and the virial expansion developed in 
earlier sections. In section VI, we conisder a hard disc 
system constrained to live on a negatively curved triply 
periodic surface such as those found in the cell or among 
diblock copolymers. We summarize our results in section 
VII. 



II. TILINGS AND DISC PACKINGS 

We begin our study of the hyperbolic plane, and curved 
space in general, with a discussion of the classical prob- 
lem of tiling a plane with regular patterns. Regarding 
the hyperbolic plane in particular, we outline the fea- 
tures of our chosen model, the Poincare disc model, in 
the Appendix. 

The problem of understanding the possible coverings of 
a fiat plane with uniform tiles is deceptive in its simplic- 
ity, and the different shapes the tiles may take or patterns 
they may be adorned with have long been pondered. In 
the last 100 years or so it has become known that there 
are precisely 17 ways of covering a flat plane with these 
tiles, corresponding to the 17 so-called Wallpaper Sym- 
metry Groups [131 . The Alhambra, in Granada, Spain, 
for example, is famous for its artistic representations of 
all 17 possible tihngs 0. 



A. Tesselating Two-Dimensional Spaces 

If we restrict ourselves to covering the plane with only 
unadorned polygons, as one would consider in a Voronoi 
tesselation, however, we reduce the number of possible 
tilings considerably. Furthermore, since crystalline order 
is a necessary starting point, we also restrict our tiles to 
regular polygons. We then begin with a regular lattice of 
Voronoi cells and consider a tesselation of F p-gons, with 
q meeting at each vertex. We allow for the possibility 
that F may be either finite, corresponding to a compact 



4 



object, or infinite, corresponding to the Euclidean or hy- 
perbolic plane. Each p-gon has p vertices, but each vertex 
is shared by q equivalent p-gons. Thus the total number 
of vertices is pF/q. Likewise, each edge is shared by two 
p-gons and so the total number of edges is pF/2. The 
total number of faces is F. The Eulcr character x of the 
surface constrains these three numbers and it follows that 



X = V-E + F = pF 
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(9) 



For the flat periodic plane (the torus) x = 0, so 
p~^ + q~^ — 1/2 and we flnd the three tilings repre- 
sented by the Schlafli symbol {p,q} = {6,3}, {3,6}, 
and {4,4}, hexagons, triangles, and squares, respec- 
tively. The spherical topology is the only one with pos- 
itive x{— 2), admitting the five Platonic solids {p,q} — 
{3, 3}, {3, 4}, {3, 5}, {4, 3}, and {5,3}: the tetrahedron, 
octahedron, icosahedron, cube, and dodecahedron, re- 
spectively. Finally, we turn to the pertinent geometries 
with negative Gaussian curvature, for which x < a-nd 
p~^ + q~^ < i. In this case there are an infinite number 
of integral pairs {p, q} and an infinite number of corre- 
sponding regular tesselations, for example, {4, 5} (Figure 
[T|). As in the tilings of the Euclidean plane, the {p, q} tes- 
selation is dual to the {q,p} tesselation - the vertices of 
one become the face centers of the other. However, this 
does not mean that their packing properties are the same, 
just as the packing of discs on the centers of hexagons is 
much more efficient than packing them on the centers of 
triangles. 



B. Building Ordered Disc Packings 

We use these tesselations to produce close-packed, 
crystalline-ordered arrangements of discs simply by in- 
scribing circles into their constituent regular polygons. 
Note that this procedure results in the familiar hexago- 
nal crystal when applied to the {6, 3} tesselation of flat 
space. In general, those packings produced from tessela- 
tions with q = 3 will result in the most efficient, densest 
crystals, as q — 3 corresponds to local triangular close- 
packing around each vertex. The rigidity or stability of a 
packing is another feature of interest that may be inferred 
from the parameters of the packing. In particular, a con- 
figuration is considered to be isostatic [26j if there are ex- 
actly enough contacts to balance the degrees of freedom 
in the system, that is, the system is at the threshold of 
stability. In an A'^-particle packing in d dimensions there 
are dN degrees of freedom and ^zN contacts, where z 
is the average number of contacts per particle, so the 
condition for isostaticity is z = 2d. Therefore, for an 
arrangement to be stable we must have p > 4, with iso- 
staticity at p = 4. Note that isostaticity is a topological 
issue and is independent of the curvature of the surface. 

A central property of these crystalline arrangements 
is their area fraction, as it indicates which packings are 




FIG. 1: (Color online). A {4,5} tiling of a segment of the 
hyperbolic plane, represented in the Poincare disc model, with 
inscribed discs. Despite the apparent differences, the discs 
are of uniform size and the tiles are all equivalent regular 
polygons. The red triangle displayed in the central cell is 
used to calculate (h. 



favored entropically near a given curvature. We find the 
area fraction of the {p, ^{-packing through application of 
the Gauss Bonnet Theorem and hyperbolic trigonometry. 
When there are F identical polygons, the area of each can 
be found from the Eulcr character x- 
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The radius of the incircle can be determined via the 
dual law of hyperbohc cosines [l^l, cos\i{-\/ —Kr) = 
cos(7r/g)/ sin(7r/p) (Figure [T]). The area fraction for the 
{p, (/{-packing is thus: 



cos(7r/g) 
sin(7r/p) 



1 
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(11) 



We recover the appropriate packing fractions, 7r/3\/3, 
7r/4, and 7r/2\/3, for the Euclidean tilings, {3,6}, {4,4}, 
and {6, 3}, respectively. As p — s- 00 for close-packed q = 3 
tesselations, we find (j) = 0.9549, the known packing frac- 
tion of the best packing on the hyperbolic plane at any 
curvature [1^. The lowest packing fraction for an iso- 
static configuration (p = 4) is ^{4,00} = V^— 1 ~ 0.4142, 
an area fraction far below that for the Euclidean square 
lattice 7r/4 « 0.785. Also of note is the broadening of 
the range in that supports stable configurations, from 
^ 0.13 in two-dimensional flat space to Acj) 0.54 
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III. FREE-AREA THEORY IN A UNIFORMLY 
CURVED BACKGROUND 

Now that we have exphcitly constructed close-packed 
disc configurations with crystalhne order, we have a 
foothold on the analysis of a thermodynamic fluid of hard 
discs li ving in a curved space. Kirkwood's free-volume 
theory [23| is highly successful at modeling the equation 
of state for traditional crystals, particularly near their 
close-packed, maximum density 'sO, fsTj . Somewhat sur- 
prisingly, free-volume theory remains successful at lower 
densities as well; as such, it provides a useful tool for 
probing a thermodynamic system, so long as the high 
density, high correlation behavior is known. 

Before we apply free- volume theory to our curved sys- 
tem, recall that in flat space the shapes of the Voronoi 
cell and the free-volume cell are equivalent and, thus, 
only the scaling of the free-volume cell is important to 
the equation of state. This scaling is dependent only on 
dimension, and we find a simple expression for the equa- 
tion of state in terms of the maximum packing fraction 
^max and the hypervolume per particle, v |32l |: 




K = -16 



FIG. 2: (Color online). The equation of state for free volume 
theory in a disc shaped Voronoi cell. The top, green curve 
is for flat space; subsequent curves are at increasing negative 
curvature. Warmer colors correspond to more curvature. 



Notice that we recover the expected Euclidean result 
of 7rr^(l/(/)^/^ — 1)^ as we take K ~* 0. Furthermore, 
denoting the proper area of the disc as a, we have 



Pfv = - 



ln{[((?!)max/0) 



1/D 



If} (12) 



Unfortunately, an analogously simple treatment does not 
exist in the presence of a curved background. The curva- 
ture itself, inducing geodesies to converge (positive K) or 
diverge (negative K), breaks the scale invariance enjoyed 
by all shapes on R", and hence destroys the similarity re- 
lation between the Voronoi cell and the free-volume cell. 



A. An Idealized Circular Free Area Cell 

As an alternative, we consider an idealized case in 
which the symmetry of the system itself ensures that 
similarity is restored. If, for example, we have a two- 
dimensional fluid composed of hard discs, and we posit 
that the free area cells are also circular, then we may 
proceed in direct analogy with Euclidean space. The free 
area is simply the area of a disc whose radius is the dif- 
ference between that of the Voronoi cell and that of the 
particle, Af^ec = -^{cosh[y/-K{rccn - ^disc)] - !)• The 
angle addition formula for the hyperbolic trigonometric 
functions, along with the recognition that (j) is the ratio 
of ^disc to ^cciii allows us to write: 



A, 



-[co^\\J —Kr\ + (/) — 1) cosh[\/— Xr] 



1 



-V(cosh[\/^r] - l)(cosh[V^r] + 2(/>- 1) 
X sinh[V^r] - 1 y (13) 
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(14) 



and it follows that the pressure diverges as a simple pole 
as (/) 1. Notice, however, that as we increase the 
background curvature the divergence at = 1 becomes 
sharper, apparent at ever higher area fractions (see Fig- 
ure[2]) , indicating a lower system pressure than is seen in a 
corresponding flat space system. This suggests that both 
particle collisions and more complicated, many-particle 
arrangements become increasingly scarce as the curva- 
ture rises. 



B. Free Area For Tesselation-Derived Packings 

What if we were to consider a non-idealized system 
based on the ordered packings discussed in section II, 
with polygonal Voronoi cells? In order to cope with 
the shape distortions caused by the background cur- 
vature, we divide the Voronoi cell into individual re- 
gions, as depicted in Figure [3l and calculate the ar- 
eas of the strip and corner regions separately, implying 
Afreo = Avoronoi-P^strip-Mcorner- Without loss of gen- 
erality we may set if = —1 since it is only the product of 
the disc size and \f--K which enter into all our expres- 
sions. We can access different curvatures by changing the 
size of the discs. Direct application of the Gauss-Bonnet 
Theorem yields ^voronoi = p7r — 27r — 2p7r/<?. In order 
to find the other two areas, we first need s, J, and a. 
By considering the triangle formed by the Voronoi cell's 
center and two adjacent vertices we may apply the dual 
hyperbolic law of cosines (Appendix A. 2) to find: 
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FIG. 3: Partition of a {p, q} Voronoi cell into its free area cell, 
p corners, and p strips. Note that the boundaries of the free 
area cell are not geodesies. 



{4,10«} {4,10} {4,5} {8,3} {12,3} 




FIG. 4: (Color online). The equation of state for free area 
theory for a range of {p, q} tesselations. Note that the pres- 
sure curves diverge at their maximum packing fractions and 
that they differ even at low volume fractions, 0. The vertical 
dashed lines indicate the maximum volume fraction for the 
tesselations listed on top of the figure. The (light blue) Xs 
are the result of our molecular dynamics simulation for curva- 
ture just below the {4, 5} tesselation and the (black) crosses 
are for curvature just above {4,5}. The number of particles 
ranges from 1 to 9. 
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To proceed we construct two identical triangles out of the 
corner 4-gon by bisecting a. The dual hyperbolic law of 
cosines applied to the side of length r then gives: 



a = sm 



cos - 

"3 

cosh(r) 



(16) 



Armed with an expression for a, we apply the hyperbolic 
law of cosines again, this time to the side of length S 
yields: 



cosh ^ 



cosh^(r) — cos^ 
cosh(r) sin 
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Now Astrip must be the area within a given distance 
(namely, r) of a geodesic of length s — 26. This area 
is (s — 25) sinh(r) (see Appendix A. 3). Given a, direct 
application of the Gauss-Bonnet Theorem gives ^corner- 
Finally, we may eliminate r through (/>: 



= cosh 



- — -(/)+ 1 



(18) 



Collecting terms and setting t = 1 — ^ and Tpq = 1 — 
cos(7r/g)/ sin(7r/p), we have 
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Afrcc = 2psin ^ 
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qj T 
— 2p v^T^ — 1 cosh~ ^ 



2tt 



cos(|) 
sm(f) 



(19) 



+2p\/t^ ~ Icosh ^ 



l-^sin^(^) 
sin(f) 



As a check of the complex algebra, note that as — > 
0{p.g}i ~^ '''{p-q} ^^'^ Afroe — > 0. Likewise, as <^ ^ 0, 
T — > 1 and 



''^-2--q 



2tT = -2tTPx/F = Avoronoi (20) 



Unfortunately, away from these limits, this expression is 
difficult to analyze, so we instead turn to FigurelH where 
we plot the free-area equations of state for various {p, q} 
tesselations using Equation (fT4|) 

It is important to note that by fixing the lattice of 
Voronoi cells we are forced to change the packing frac- 
tion by changing the disc size. Since this is equivalent to 
changing the curvature, our free-area equation of state 
samples different curvatures at different area fractions (j). 
If we wish to apply this free-area theory to a system of 
constant disc size to sample a constant background cur- 
vature, we must allow the Voronoi scaffolding to vary 
continuously. We will have to consider tesselations with 
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FIG. 5: (Color online). Fixed radius equations of state with 
a sampling of {p, q} tesselations as their divergent endpoints. 
The solid blue curves axe for fixed q, and the dashed red 
curves are for fixed p. Note that in almost all cases these 
curves lie directly on top of one another. The inset displays 
the parameter space for the packing parameters {p,q}. The 
thick red line is the allowed combinations of p and q in flat 
space and the purple line beneath it a representative example 
from positively curved space. The blue curve gives the allowed 
combinations for with r — 0.53. The {4, 5} packing is 
identified and the paths through parameter space taken by 
holding either p or g constant are shown. 



non-integral p and q\ we can interpret these as the aver- 
age number of edges per polygon and the average number 
of polygons meeting at a corner in a tesselation made of 
a variety of polygons. This is in the spirit of recent work 
on average polyhedra [33*1 which have a non-integral num- 
ber of sides. To maintain the curvature as we vary (p, we 
might fix p or q adjusting the other so as to keep r fixed in 
(fT5)l . We might also hold p/q fixed and set p to maintain 
r. The freedom afforded by choosing one's path through 
p-q parameter space is troubling - there seems to be no 
a priori way for one such path to be chosen over another, 
and we are forced to wonder how the behavior of the sys- 
tem is determined. Fortunately, upon considering specific 
cases, the problem disappears, as the fixed-radius equa- 
tion of state displays insensitivity to this choice (Figure 
[5]) , particularly near high density, indicating that perhaps 
there is a hidden, unaccounted-for symmetry among the 
isotropic Voronoi scaffoldings. We only find path depen- 
dence for very large values of q and p - there are no real 
solutions for p when q is large until we get close to the 
divergence in the pressure, though there the divergence 
is insensitive to the details again. 



IV. THE VIRIAL EXPANSION TO THE 
EQUATION OF STATE 



On the other hand, we are also interested in the system's 
behavior when it is much more sparse. The expansion 
about low densities provided by the virial expansion to 
the equation of state 3J| is suitable in this regime, al- 
though the connection between the two regimes is still 
not well understood in hard-core systems [35|. In con- 
trast to free-area theory, the virial expansion's validity 
is dependent on the ability to take a proper thermody- 
namic limit. We therefore resort to the periodic spaces 
discussed earlier as our solution to this difficulty. 

We note, however, that for hard discs, the virial ex- 
pansion only includes local interactions among clusters 
of particles. In fact, for pure, hard-core potentials, the 
successive terms of the virial expansion have a simple in- 
tuitive interpretation. The first order term, as always, 
is simply the ideal gas law. The next term removes the 
contribution from particle pair configurations that are 
disallowed by the steric interaction, that is, pairs of par- 
ticles whose hard-cores overlap. But this over-counts, as 
the same configuration with three or more particles si- 
multaneously overlapping is subtracted multiple times. 
The third virial term explicitly corrects this problem for 
three-particle arrangements, but similarly over-counts for 
larger numbers of particles, this time in the other direc- 
tion. Higher terms continue to both correct over-counting 
errors from previous terms and contribute smaller over- 
counting errors of their own, in a manner reminiscent 
of an infinite tower of Venn diagrams. Fortunately, as 
long as the local clusters under consideration do not wrap 
around the periodic space and self-interact, we may ig- 
nore the topology introduced by our usage of a periodic 
space and perform the cluster integrals on H^. In or- 
der to ensure that this is possible, we may continue to 
increase the area available to a cluster by constructing 
ever higher genus manifolds [22]. Furthermore, for man- 
ifolds with constant K <0, Zi^ A[\'j,'^ + K/{12n)] > 
when < —\2'it/K. In other words, there is an ample 
range of K in which to work where the curvature of the 
manifold does not probe the quantum regime of the gas. 
Hence, the pressure P = NksTdhi Z / dA is independent 
of the proportionality constant. 

Just as in the last section, we take advantage of the 
lack of scale invariance on to probe higher curvatures 
not by changing the ambient space, but simply by consid- 
ering larger and larger discs. We consider discs of radius 
r\J —K. Without loss of generality we again fix the Gaus- 
sian curvature to be X = —1. Dimensional analysis may 
always be employed to insert the appropriate factors of 
K. As in fiat space, the second Virial coefficient, B^ir), 
is given by half the excluded area of a single disc. To see 
that this is true, recall that, in general, Bi(r) is given by 



B^ir) = — - 



fi2 d^ri d^ra 



(21) 



Thus far, by considering the extension of free- area the- 
ory to uniformly curved spaces we have been working in 
the high-density regime, since free-area theory is the first 
order term in an expansion about a close-packed state. 



where /12 is the Mayer function and the position integrals 
for both particle 1 and particle 2 are carried out over all 
space. Since is isotropic and homogeneous, we may 
integrate out one particle position and pick up a factor 
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of V. The remaining integration is over the separation 
between the discs. Meanwhile, the Mayer function for 
hard discs in the Poincare disc model is: 



— 1 if cosh 
otherwise 



1 



< 2r 



(22) 



so long as one of the discs is at the origin. The expression 
when both discs are away from the origin is more compli- 
cated but follows from (|A4|1 . However, for the purposes of 
B2(r), we can always put one disc at the origin. Turning 
the step function into a condition on the bounds of the in- 
tegral over X, inserting the metric from the Poincare disc 
model (see Appendix A) , and integrating out the angular 
part results in the following expression for B2{r): 



dx 



(1 



(23) 



Notice that this is precisely the circularly-symmetric, hy- 
perbolic, area integral needed to calculate the area of a 
disc of radius 2r on H^, and indeed, carrying out the in- 
tegration leaves us with B2{r) ~ n [cosh(2r) — 1], which 
is half the excluded area of a single disc, the classic result 
from the virial expansion. Further note that as r — > 0, 
B2 ~ 27rr^, we recover the d = 2 Euclidean result; in d 
dimensions B2{r) = (Att)^ r'^ /[r{f)d]. 

We may carry out a similar procedure to calculate 
higher coefficients. The general relation for Bz{r) is: 



Bz{r) 



1 



/12/23/13 d'ri d^ra d^v^ (24) 



and we may still integrate out one of the disc positions. 
Unfortunately, regardless of which disc we set to the ori- 
gin, we are now guaranteed to retain a Mayer function 
where neither of the discs whose interaction it represents 
have had their position integrated out. This means that 
we are forced to contend with an integral of the form 



B. 



327r r, 

Xi 



2ti I'b 

d9i2 / dx2 
Jo 



(25) 



X2 



{1-xir (1-4) 



2h2{xi,X2,9i2) 



where h is the same upper limit as in Eq. 1231 and the 
remaining Mayer function is written in terms of the co- 
ordinate distance of the two particles from the origin and 
the angle between them: 



-1 



J12 



1 



2^ 



2xiX2 COS 612 



otherwise. 



(1 



1 



< cosh(2r) 



In order to evaluate this integral, and the similar ones 
that appear in the higher virial coefficients, we turn to 
numerical, Monte Carlo techniques. In addition, as is 
traditional in flat space ^36,], the higher coefficients -B„ 
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FIG. 6: (Color online). The second, third, fourth and fifth 
virial coefficients in appropriate units of the area, a, as func- 
tions of disc size (for K = —1) r. The zero size limit matches 
the fiat-space coefficients. B2 is exact and error bars are in- 
dicated for all other coefficients. These errors are given by 
the estimated variance of the mean of the integrand, and, for 
large 77, scale like 77"^''^, where 77 is the number of Monte Carlo 
trials. 



are reported in units of a"~^ — [A(r)]" to suppress 
size effects and generate the expansion for P/{pkBT) in 
the area fraction, <j) = pa: 



P = pkeT 



E 



(26) 



We have checked the r limit in all cases, and find 
agreement with the known Euclidean results. We plot 
B2/a, B^/a?, B^/a^, and B^/a^ in FigureElas functions 
of r. To the accuracy of our calculation, over the range 
of curvatures considered, iJa, i34, and B^ all remain pos- 
itive. There is evidence that, in higher dimensional fiat 
spaces, the virial coefficients alternate in sign [s^ [s^l 
and, hence, the leading singularity that controls the ra- 
dius of convergence of the expansion sits on the negative 
real axis [s^l • The location of this singularity, and thus 
the range of applicability of the virial expansion itself, 
remains an open question in D = 2, 3. Though not con- 
clusive, our results here suggest that the introduction of 
curvature does not qualitatively change the behavior or 
location of the leading singularity in the expansion. 



A. Pade Approximants and the Radius of 
Convergence 

The quantitative behavior of the virial expansion's di- 
vergence as a function of the curvature, on the other 
hand, may be extracted from the coefficients that we 
have calculated, at least as a trend. To this end we have 
approximated the virial series by ratios of rational func- 
tions, the Fade approximants [l^. We find the coeffi- 
cients of two polynomials am{p) and bn{p), of degree m 
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and n, respectively so that 

p^..rJ. + gB,yj^:|^ (27, 

The left hand side of (P7|) is fifth order and so there are, 
in principle, six arbitrary coefhcients. As a result, we can 
only study the Pade approximants with n + m < 5. We 
look for the divergences in these approximants to esti- 
mate where the virial expansion diverges - an estimate 
of the ultimate packing fraction if the virial expansion 
converges for all physical p. Our results are shown in Ta- 
ble [D Notice that, for the Pade approximants with only 
one pole (that is, of the form [to, 1]) there is a clear trend 
towards divergence at lower p as the curvature increases. 
Furthermore, it is worth pointing out that among these 
approximants, the one that is the most sensitive, the 
[3, 1], also agrees best with the putative location of the 
divergence of the fluid branch in flat space, at around 
(p ~ 0.84, as r goes to zero. 



[m, n] 




r ^ 


r = 0.5 


r = 1.0 


r = 1.5 


r = 2.0 


[0,4] 


-3.444 


4.335 








[1.3] 


2.720 


-0.592 


-8.695 


0.797 


0.650 


[2,2] 








0.798 


0.695 


[3,1] 


0.798 


0.797 


0.794 


0.770 


0.676 


[0,3] 


-5.145 


-45.09 


0.634 


0.409 


0.258 


[1,2] 






0.888 


0.821 


0.829 


[2,1] 


0.735 


0.733 


0.732 


0.730 


0.737 


[0,2] 


0.737 










[1,1] 


0.639 


0.637 


0.630 


0.620 


0.611 


[0,1] 


0.500 


0.470 


0.393 


0.298 


0.210 



TABLE I: The leading divergences on the real axis in </!> at var- 
ious values of curvature and for different Fade approximants 
about the origin. Empty entries indicate that all poles for 
that approximant are located off the real axis. Though (p has 
physical meaning only in the range [0, 1], a divergence at an 
unphysical value of (j> may still be meaningful, if, for example, 
it encodes information about the radius of convergence of the 
virial expansion. 



V. MOLECULAR DYNAMICS ON THE 
HYPERBOLIC PLANE 

Equipped with these theoretical models for the system, 
particularly in the high and low density regimes, we turn 
to the full-blown hard-disc fluid. We choose to employ 
constant energy molecular dynamics. There is some sub- 
tlety in handling the boundaries of the simulation region: 
recall that the area near the edge of a shape scales at the 
same rate of the area of the entire shape on H^. This 
means that the effect of hard boundaries cannot be min- 
imized by considering larger regions. We thus choose to 
implement periodic boundary conditions, although there 
is subtlety here as well. As mentioned above, because 



of the scale dependence of the hyperbolic plane and the 
topological requirements of the periodic region, we are 
restricted to choosing for our simulation region the unit 
cell from tesselations of the form {An, 4n} for n an in- 
teger greater than 1. This ensures both that the sides 
of the polygon may be identified in a systematic way to 
produce an orientable manifold and that, upon identifi- 
cation, there is not a conical singularity at the vertices 
of the polygon (so there are exactly 2tt radians of space 
around each vertex). 

We have chosen not to use the event-driven molecu- 
lar dynamics traditional for pure hard core systems, as 
projecting the next collision is extremely time intensive 
when the equations of motion are nontrivial as they are 
in a curved background. Instead, we evolve the system 
by time steps, checking for and resolving collisions as 
they occur. The collisions themselves require some care 
as well. The center of mass frame in a curved space 
is not necessarily inertial, even if the constituent parti- 
cles are all following geodesies - imagine a sphere with a 
particle moving around the equator and another sitting 
motionless on the North pole - the center of mass follows 
the 45^^ parallel which is not a geodesic. As a result, 
the textbook techniques and formulae for handling colli- 
sions do not apply. As an alternative, we decompose each 
particle's velocity into the component along the geodesic 
connecting the particle centers at contact and the compo- 
nent normal to this geodesic. Note that these directions 
will differ between the two colliding particles in the lo- 
cal coordinate system due to the complexity of parallel 
transport on H^. In order to resolve the collision, the 
components of velocity along the contact geodesic are 
swapped, and the other components remain unchanged. 
Hence, collisions are guaranteed to conserve energy and 
momentum. 

To calculate the equation of state from our simulation, 
we employ the collisional virial expression for the pressure 

P = ^(NkT+^ ^ ^I^pI) 

\ collisions / 

where A is the area of the simulation region, r is the time 
over which collisions are summed, d is the disc diameter 
and Ap is the momentum transferred from one particle 
to another during the collision. Finally, we remind the 
reader that to fill out the equation of state with differ- 
ent values of cf) it is necessary to change the number of 
particles simulated, TV, as varying the disc size means 
changing the effective curvature. 

A. Systems Near the {4, 5} Tesselation 

We probe curvatures r = 0.53 and r = 0.531, just 
below and just above that for the {4, 5} tesselation, 
^{4,5} = 0.5306. At these curvatures no regular tessela- 
tion is allowed, and so non-crystalline arrangements give 
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the best packings [sP]. Thus we expect to probe dis- 
ordered configurations even at high densities with this 
choice. At lower densities, however, the gas will still be 
sensitive to the {4, 5} tesselation and should follow the 
free-area theory prediction for the {4,5}, failing only at 
higher area fractions. In particular, we expect that the 
system's pressure should be lower than that of the iso- 
static crystal for higher and diverge at larger 0. This is 
precisely the case in three-dimensional flat space where 
the liquid branch is lower in pressure than the the free- 
volume crystalline pressure at low volume fractions. In- 
deed, these expectations are borne out by the data (Fig- 
ure H]). 



B. Approaching Flat Space from Below: 
Monodispersity and Disorder 

In addition to investigating the system at packing frac- 
tions close to an isostatic crystal, we might use a limiting 
procedure to take the background curvature to zero. By 
looking at the fluid in low curvatures, we expect to gain 
some insight into hard discs in flat space. In particu- 
lar, looking at curvatures near, but not at, zero guaran- 
tees that we will avoid crystallization, as no crystal pack- 
ings are available in this curvature range. It is therefore 
possible that we may generate disordered, monodisperse 
conflgurations that can be carried into flat space with 
low probability of introducing overlaps or defects. Un- 
fortunately, due to the broken scale invariance on H^, a 
very large number of simulated discs is needed to probe 
this regime, making the simulation numerically unten- 
able, even for the high-speed digital computers of our 
day. 



VI. PACKING AND DYNAMICS ON 
NEGATIVELY CURVED SURFACES IN FLAT 
SPACE 

Thus far, we have used uniform, negative, background 
curvature as a way of constructing an idealized model 
in two dimensions for a more complicated three dimen- 
sional system. Owing to the fact that no surface with 
uniform, negative curvature can be isometrically embed- 
ded into R^, however, this model is unphysical. But 
are there examples of real systems that are qualitatively 
close, for which the predictions of our model may be 
brought to bear directly? Triply periodic minimal sur- 
faces, with their numerous physically advantageous prop- 
erties [i^l, their prevalence in smectic liquid crystals, di- 
block copolymers, and even biology, and their status as 
inflnite surfaces with negative curvature everywhere are 
ideal for analysis with insights from the abstract hyper- 
bolic plane. 

In lieu of a direct molecular dynamics simulation, as 
done on H^, we instead take advantage of the numeri- 
cal insight provided by examining allowed and disallowed 



conflgurations [4l|. In fact, we may obtain the pressure 
simply by randomly sampling particle conflgurations at 
a given density and determining the fraction that violate 
the no-overlap condition of the hard-core potential. If 
the probability of not flnding an allowed state at density 
<j) -f (50, given an allowed state at 0, is J{(t>)6(l), then the 
following relation holds [4]| : 

Ji'l')^Tr^m)-'I^T/v] (29) 

where v is the area per disc, V the system volume, T the 
temperature. 

Hence, this approach gives us the equation of state. 
We have traded the numerical difficulty of dynamics of 
particles constrained to a complicated surface with the 
algorithmic difficulty of finding a way to generate points 
at random on an arbitrary surface such that the distri- 
bution is flat with respect to the local area. 



A. Randomly Distributed Points on an Arbitrary 
Surface 

We have employed a novel geometric algorithm for se- 
lecting points from an arbitrary surface randomly with 
even distribution in the area (so long as it is either com- 
pact or periodic in each dimension). Knowledge of the 
metric is not required, only a set of parametric equations 
that deflne the surface in M."^. The algorithm works as fol- 
lows. First, one repeat unit of the surface in whichever 
directions are periodic, and the entire surface in those 
directions for which the surface is compact, are consid- 
ered to be inside a box. At this stage it is necessary 
to determine the maximum number of times a straight 
line can intersect this region of the surface, say /max ~ 
this can typically be done analytically, though for a very 
complicated surface numerics may prove more efflcient. 
The box is inscribed in a sphere, and a great disc is cho- 
sen at random in this sphere. Next, a point is chosen 
at random on this great disc and the line normal to the 
great disc is drawn through it. Finally, we determine the 
intersections of this line with the surface and randomly 
select one from among them. To ensure that the result- 
ing distribution of points is flat with respect to area, it 
is necessary to begin the entire process anew if the line 
fails to intersect the surface. Furthermore, if the num- 
ber of intersections is less than /max then an appropriate 
number of placeholders must be included when randomly 
selecting an intersection point. If a placeholder is selected 
then the entire procedure is repeated. If an actual inter- 
section point is selected, then that point is the output of 
the pass through the algorithm. 

To prove that this procedure yields a distribution flat 
with respect to the local area, consider the shadow cast 
by the surface when light from an arbitrary direction 
strikes it. Now consider the contribution to this shadow 
from a differential piece of the surface. The area of the 



11 



part of the shadow cast by this small patch will be di- 
rectly proportional to the likelihood of selecting a point 
from it when the great disc chosen at random is normal 
to the direction of light, P{dA) oc dAj cos((?)|, where 9 
is the angle between the patch's surface normal and the 
light rays. But the algorithm specifies that we constantly 
choose, at random, new directions for the illumination. 
This has the effect of averaging over all incoming direc- 
tions for the light, 9 and (j), and removes the dependence 
on the angle, leaving the likelihood of selecting a point 
from any given patch of surface directly proportional to 
only the area of the patch. To see this, imagine a par- 
ticular great disc parameterized by x and y. The surface 
has a shadow on the disc of area Ashadow. Ignoring, for 
the moment, multiple leaves we may parameterize the 
surface in terms of x and y so that the probability of 
selecting any particular area element dA given that we 
pick a point in the shadow is: 

P{dA I the point is in the shadow) = (30) 
_ dA\ cos 6*1 _ dA\ cos 6*1 

/shadow '^^^y ^ 

shadow 

Note that ^shadow depends on the choice of great disc, so 
this conditional probability does not pick area elements 
with equal weight. However, the probability of picking a 
point in the shadow is just AshadowMgreat disc- Putting 
this together we see that P{dA) = dA\ cos0|/Agroat disc- 
Integrating over all directions yields a flat distribution. 
Finally, if there are multiple leaves of the surface behind 
one another with respect to a given direction of incident 
light, randomly selecting from among them and multi- 
ple placeholders, as described above, guarantees that no 
patch is penalized or rewarded for its position relative to 
distant pieces of the surface. 

B. The Plumber's Nightmare 

The simplest triply periodic minimal surface, and the 
one that we consider first, is the Schwarz P-Surface, 
sometimes referred to as "The Plumber's Nightmare" for 
the impression it gives of a never-ending maze of pipe 
work. This surface appears in the study of lipid bi-layers 
and liquid crystals [i^. In lieu of using the ex- 
act minimal surface, however, we use the trigonometric 
approximation, as its departure from minimality affects 
only the mean curvature, leaving the Gaussian curvature 
nearly unchanged. Our approximate Schwarz P-Surface 
is defined by the following relation: 

cos(x) -I- cos(y) -I- cos(z) = 0. (31) 

Using this level set approximation, we have used our 
algorithm to choose points as depicted in Figure [71 We 
have checked our method by projecting these points along 
an axis and normalizing by the square root of the metric 
to see a uniform density. With this in hand we may now 
calculate the pressure using ([M)) . We lay down N discs of 
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FIG. 7: (Color online). Dimensionless pressure as a function 
of A*' for discs of various sizes living in the P-Surface. The disc 
diameters range from 0.05 to 0.5 in units where the repeat 
length of the P-Surface cell is 2tt. Warmer colors correspond 
to smaller discs. A subset of the random points used is shown 
in the upper left corner. 



a given size, one at a time, checking to see if they can be 
inserted without overlapping an existing disc. To accept 
or reject each disc, we check the distance between points 
along the surface and compare this to the disc diameter. 
The distances between disc centers along the surface are 
also estimated in a metric-free way, by taking a small 
step along the line in M'^ connecting the disc centers and 
finding the plane normal to this line at the location of the 
step. We then calculate the curve of intersection between 
the P-Surface and the plane and find the point on this 
curve closest to the initial disc center. Tabulating this 
distance and then repeating the process with the location 
of the last step as the new initial point from which the 
next line in R"^ is drawn completes the algorithm. Armed 
with an unambiguous state acceptance criterion, we have 
an avenue to calculate J{4>), however, our relevant control 
parameter is N. So we instead write: 

J{N) - ^ [P{N) - NT/V] (32) 

using (p — Nv/V. Isolating a dimensionless pressure, 
P = P[N)V/T allows us to wash out the hard core sys- 
tem's insensitivity to the temperature and the overall 
scale invariance, leaving: 

P = n(i + ^J{N)). (33) 

The dimensionless pressure calculated in this way is 
shown in Figure [7] as a function of N . Notice that for 
all considered disc sizes, the pressure shows little devi- 
ation from an ideal gas state {P — N) until near the 
divergence. This stands in excellent agreement with our 
predictions from free-area theory and meshes well with 
the observations made on hard core systems in general 
by Kamien and Liu [4l| . 




12 




100 200 300 400 500 600 

N 

FIG. 8: (Color online). Dimensionless pressure as a function 
of A*' for discs of various sizes living in the D-Surface. The 
disc diameters range from 0.05 to 0.5 in units where the repeat 
length of the D-Surface cell is 27r. Warmer colors correspond 
to smaller discs. A subset of the random points used is shown 
in the upper left corner. 

C. Chaos Carolinensis and the Schwarz D-Surface 

Another triply periodic minimal surface of interest is 
the Schwarz D-Surface, or Diamond Surface. Lipid bi- 
layer and liquid-crystalline systems also both exhibit this 
surface, though it appears in a more surprising context 
as well - biology. When the amoeba Chaos Carolinensis 
is deprived of food, its mitochondrial membrane tran- 
sitions from an amorphous, unpatterned state to a D- 
Surface, possibly to optimize its surface-to-volume ratio 
while cannibalizing its own membrane phospholipids for 
energy [2^. We again choose to use the trigonometric 
approximation; this time the defining relation is given 
by: 

cos{x — y) cos(z) -I- sin(a:; + y) sin(z) = 0. (34) 

As before, we generate points on this surface (see Fig- 
ure [S]), and again we calculate the dimensionless pressure. 
For the full range of considered disc sizes there is little 
deviation from the ideal gas until near the divergence, 
similar to the P-Surface, as shown in Figured Perhaps 
this is not surprising, as the fundamental cell of the P- 
and D-Surfaces are made of the same local pieces patched 
together in different ways. Thus the purely local interac- 
tions are likely insensitive to this longer range structure 
leaving open the possibility that only interactions on the 
scale of the fundamental cells might be sensitive to the 
changed topology. 

VII. DISCUSSION 



brings the qualitative features of the system more in- 
line with the frustrated three-dimensional flat case than 
with the un-frustrated two-dimensional system. We pre- 
sented theoretical models of the system both at low and 
high density, with the Virial expansion and the enumer- 
ation of crystalline states along with a free-area theory 
based around these crystal endpoints. We simulated the 
system in a periodic piece of , and compared these re- 
sults with our theory and with the well-known flat-space 
systems. Finally we numerically analyzed a similar sys- 
tem on a physically realizable curved surface with non- 
constant negative curvature. 

We emphasize, however, that these are only the first 
steps into a rich geometric framework. Technological ad- 
vances will allow a clean look at the system as the cur- 
vature is turned off and will make higher Virial coeffi- 
cients attainable. More complicated particle shapes and 
interactions will allow a complete modeling of sub- micron 
scale packings constrained to curved membranes, as ex- 
ist in lipids, diblock copolymers, and abound in biology. 
In addition, it is worth considering and the delicate 
balancing of factors that allow the thermodynamic limit 
to be trivially sensible in flat space and difficult to grasp 
anywhere else. 
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APPENDIX A: CLASSICAL NON-EUCLIDEAN 
GEOMETRY 

Here we review some properties of and the well- 
known Poincare model, and collect some useful classical 
identities and theorems on uniform curved geometry in 
general 

1. The Poincare Disk Model: Metric and Distance 

We employ the Poincare disc model of H-^, popularized 
by Escher 44], which conformally maps the hyperbolic 
plane with curvature K to the disc of radius y/—K^^ on 
R^. Using the coordinates {x,y), the metric is |22| 



In summary, we have constructed a monodisperse, two- ^^2 _ 4 (^dx + dy ^ (Al) 
dimensional system with hard-core repulsions by intro- [1 -|- if (x^ -I- y^)]^ 
ducing negative background curvature to frustrate lo- 
cal close-packing. This suppresses crystallization and wheie x^ + < —1/K. 
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Geodesies in this model are arcs of circles which in- 
tersect the bounding circle normally; circles in the hy- 
perbolic plane remain circles in this model, but their co- 
ordinate centers and radii vary, shifting outward from 
the apparent center and shrinking far from the origin, 
respectively. The radius of curvature sets a lengthscale 
and, consequently, any system on is manifestly not 
scale invariant. 

In addition, it is useful to understand that the Poincare 
disc model arises via stereographic projection [i^ in 
the same way as spherical mappings, and, in particular, 
maps of the earth. Define the following quadratic form: 



Q(x) 



-y^ + z^. Then the collection of points with 



norm one with respect to this quadratic form define a 
two-sheeted hyperboloid. We carry out the stereographic 
projection by mapping points on the sheet at positive z 
onto the unit disc in the x-y plane via the collection of 
lines passing through the point (0,0, —1). Given a point 
in the unit disc, {ux,Uy) we now have the map to the 
hyperboloid: 



2Ut 



2u„ 



1 



1 



1 



1 



(A2) 



where u is the length of the vector in the x-y plane. The 
distance function for the Poincare disc model is then 
given by the symmetric bilinear form related to Q by 
polarization: 



cosh(d(u, v)) 



i,v) = i(Q(u + v)-g(u)-Q(v)) 

(A3) 

where we have denoted the two-dimensional vectors in 
the disc by u and v and their images on the hyperboloid 
in by u and v. Substituting the expressions for the im- 
ages in terms of the model vectors and simplifying yields 
the distance function for the Poincare disc model with 
K = -1: 



1 



d(u, v) = cosh ^ 
Were we to write u = v 



u- 



ds 



2 - [d(u,v)]^ = cosh"^ 



du then we find 



(A4) 



1 



(1 -fi2)(l -52) 



(l-t22)2 

in agreement with (lAl 



(A5) 



2. The Hyperbolic Law of Cosines 

Consider a triangle in the hyperbolic plane with side 
lengths a, 5, c and corresponding angles a, /?, 7. Call the 
coordinate vertices A, B, and C. As before, we set K — 
— 1 in order to simplify notation. Now, because IHI2 is 
isotropic and homogeneous, we may assume that vertex 



A sits at the origin of the Poincare disc model. The 
distance function (|A4[) then gives the following relations: 



C = , h-^$^ and B = . h-^^^. (A6) 



cosh(6) + 1 



cosh(c) -I- 1 ' 



At the same time, applying the distance function in a 
similar way to the side of the triangle without a vertex 
on the origin gives: 



cosh(a) = 1 + 2 



IB -CI' 



(1-B2)(1_C2)- 



(A7) 



But, note that A, B, and C define a Euclidean, coordi- 
nate triangle as well. We apply the traditional law of 
cosines here: 



|B - C|2 = B^ + C^- 2BCcos{a) 



(A8) 



Substituting and ([Mjl into (K7\ eliminates B and 
C. Clearing denominators and repeated application of 
the hyperbolic trigonometric identities yields the hyper- 
bolic law of cosines: 



cosh(a) = cosh(6) cosh(c) — sinh(&) sinh(c) cos(a). (A9) 

Alternatively, there is a dual form of this relation that 
one can obtain from the cyclic permutations of the hy- 
perbolic law of cosines together with more trigonometric 
identities (hyperbolic or otherwise): 



cos(q;) = cos(/3) cos(7) -I- sin(/?) sin(7) cosh(a) (AlO) 



3. Area Within a Given Distance of a Geodesic 

Consider a quadrangular region in the hyperbolic plane 
bordered on one side by a geodesic of length d, and hav- 
ing geodesies of length r as the sides adjacent to the 
first. Finally, define the fourth side as the curve of con- 
stant geodesic curvature such that each point along it is 
a distance r from the length-d geodesic [4aj. Recall the 
geodesic deviation equation: 



(All) 



where ^(s) is the separation between two geodesies as a 
function of the arclength parameter along each, s. Con- 
sider the two sides of length r in the context of this equa- 
tion. The initial conditions are ^(0) = d and ||(0) = 
because the initial curve of length d is itself a geodesic 
and so first variations of its length vanish. We find 
^(s) = d C08h{y/ — K s) . The area enclosed by our strip 
must then be given by: 



A^ 



ds <e (s) 



(A12) 
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A = 



as required. 
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